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ABSTRACT 

Intrinsic luminosities are related to rotation velocities of disk galaxies by Tully Fisher 
(TF) relations. The Baryonic TF (BTF) relation has recently been explained with Dark 
Matter and Newtonian Gravity as well as Modified Newtonian Dynamics (MOND). 
However, recent work has pointed out that the currently used BTF relation ignores 
the contribution from hot gas and oversimplifies complex galaxy-scale physics. In this 
Letter, we advocate the use of the Neutral Hydrogen TF (HITF) relationship, which 
is free from dust obscuration and stellar evolution effects, as a clean probe of gravity 
and dynamics in the weak field regime. We incorporate the physics of hot gas from 
supernova feedback which drives the porosity of the Interstellar Medium (ISM). A sim- 
ple model that includes supernovae feedback, is generalized to include a parametrized 
effective gravitational force law. We test our model against a catalogue of galaxies, 
spanning the full range of disks from dwarf galaxies to giant spirals, to demonstrate 
that a Kennicutt-Schmidt (KS) law for star formation and simple Newtonian grav- 
ity is adequate for explaining the observed Hi scaling relations. The data rules out 
MOND-like theories, within the scope of this model. 
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^f; ■ 1 INTRODUCTION 



If all matter is accounted for, Newtonian gravity and dy- 
namics should provide an adequate description of the mo- 
tion of celestial bodies in the weak gravity regime. Yet, 
IZwickvl ([l937) found that the luminous matter in galaxies 
(then called Nebulae) could not account for the dynamics 
in galaxy clusters. This Dark Matter (DM) problem was 
extende d to galaxies with t h e advent o f reso lved rotation 
curves (|Rubin fc FordI Il970l ). iMilgromI (|l983l ) proposed a 
Modified Newtonian Dynamics (MOND), as an alternative 
to the DM paradigm. Galactic dynamics provide us a way to 
test the nature of gravity in the weak field limit and settle 
one of the most important outstanding questions in physics. 

iTullv fc Fisher! |l9z3) (TF) found a correlation between 
the global neutral hydrogen (Hi) profile widths and absolute 
optical magnitudes of disk galaxies. This relationship has 
since then, been w ell calibrated against Cepheid distances 
20001) and widely used as a distance indicator 
I2OOO'). absolute magnitudes are es- 
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sentially a mass proxy and the profile widths are a velocity 
proxy, this relation implies a connection between the mass 
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and dynamics. iMcGaugh erabl (|2000l ) have shown that even 
though gas dominated galaxies with a low stellar fraction 
deviate from the TF relation, they fall on a "more funda- 
mental" relation between the tot al baryonic mass and the 
rotation velocity. iMcGaughl l|201ll ) has shown that the BTF 
relation in gas ric h galaxies i s cons istent with the predic- 
tion from MOND. iMcGaughl (120111 ) has included only gas 
rich galaxies in the analysis, to minimize the uncertainty in 
the baryonic mass introduced by the error in stellar masses. 
However the selection by gas fraction restricts the range of 
galaxy masses and rotation velocities. This may also bias 
the slope as t he gas fraction is a function of the total mass. 
iLarsonl |20lO) has also explained the BTF, but within the 
framework of Newtonian gravity. Hence, the BTF can plau- 
sibly be explained in both Newtonian and MOND gravities, 
using a suitably selected model. Any such test is therefore 
intrinsically model dependent. 

An attempt to explain a TF relation must also take into 
account, instabilities in the disk (Dalcanton ot al. 1997), re- 
sulting s tar fo rmation and supernova feedback (Silk 199J|). 
iGnedinI |2oi3) has pointed out that BTF only accounts for 
stars, molecular, and atomic gas, while the contribution of 
ion ized gas is almost univers ally missed. Such a situation has 
led [Foreman fc Scotj (|201ll ) to comment that the claimed 
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Cold DM prediction for the BTF, bv lMcGaughl (|201ll '). is a 
gross oversimplification of the complex galaxy-scale physics 
involved. Our work, for the first time, accounts for the role of 
supernova d riven hot ga s in determining the Hi scaling rela- 
tions. In fact [sii3 |l993) has explained the TF relation in the 
B-band, which traces recent star formation, by considering 
self regulated star formation. However, the optical photons 
are strongly affected by dust obscuration and inclination, 
which may significantl y bias the TF relations. Using Near 
Infra red TF relations ( Giovanelli et al.lll997l : [Masters et al.l 
I2OO6I ) can minimize dust effects and provide tighter correla- 
tions. The alternative is to use a relation between Hi masses 
and Hi profile widths, which we henceforth refer to as the Hi 
TuUy Fisher (HITF) relation. This relation is far less sensi- 
tive to obscuration as the universe is largely trans parent at 
radio wavelengths (|Khedekar fc Chakraborti|[201ll) . The Hi 
signal originates from the ISM and is therefore independent 
of the stellar ages in the galaxy. 

Any model which expl ains this relation must 



also explain llChakrabortil 2011 ) the o bserved 



(jRosenberg fc Schneider! l2003l : iBeeum et al.ll2008bh scaling 
relation between Hi masses and areas of disk galaxie s. In 
this L etter we further develop an analytic model by ISilkl 
ljl997t ) and IChakrabortil { 201 Ih to derive the Hi masses of 
disk galaxies, taking into account the competition between 
gravitational instabilities and mechanical feedback from su- 
pernovae. We also generalize the model to include arbitrary 
(v'M /ry laws of gravity and test it against our catalogue 
of galaxies. The HITF data spans the full range of disks 
from dwarfs galaxies obser ved with the Giant Metrewave 
Radio Telescope (GMRT) tee gum et al]|2008bt). to gian t 
spirals observed with the Arecibo (iGurovich et aL 201Cf). 
The H i mass vs area data is from [Rosenberg fc Schneider! 

We conduct a joint likelihood analysis of these data 
sets within the scope of our model to demonstrate the 
validity of Newtonian gravity at galactic scales, and also 
rule out MOND. 



2 HI SCALING RELATIONS 

Hi scaling re lations in galax ies have been known for the past 
two decades. iBriggs fc Rad (1993) had already found a weak 
scahng {AV oc Mh"'^) of the velocity width (AV) with the 
Hi mass (Mhi). This however implied a strong dependence 
of Mhi on AV with a power law index between 2.1 and 
3.6 l|Briggs fc Rad 119931 ). The optical TF relation, which 
meas ures stellar mass, bre aks down at the low luminosity 
end (|McGaugh et al.ll2000l ). However no such breaks have 
been reported in the Hi scaling. We are now in a position to 
extend the HITF to the full range of galaxy masses for the 
first time and show in this work that it is consistent with a 
single scaling relation. 

Scaling relations between Mm a nd the surface area's 
of gal axies also have a long history. iGiovanelli fc Havned 
l|l983r ) reported that Hi sizes and masses were found to 
be correlated in the same manner, irrespective of whether 
they were Hi rich or Hi poor. O ptical sizes and Hi 
mass e s were also seen to be correl at ed (IHavnes fc Giovanelli 
1 1984 IVerheiien fc Sancisil I2001I ). [Rosenberg fc Schneidei 



a nearly constan t average Hi surface density of the order 
~ IO'^Mq kpc~^ [Chakraborti[ l[201ll ) has recently explained 
this surprising correlation, as the result of self-regulated star 
formation, driven by the competition between gravitational 
instabilities in a rotationally supported disk and mechanical 
feedback from supernovae. 



3 GENERALIZED MODEL 



Following [Sii3 (1993), we consider a rotationally supported 
gas disk, fragmenting and forming stars due to gravitational 
instability. In such a scenario the asymptotic velocity can be 
related to the enclosed dynamical mass (Afd). However this 
requires the knowledge of an effective force law for gravity in 
the weak field limit. Since we wish to determine the nature of 
gravity at these scales, it is necessary not to use a particular 
variant, but a parametrization which can encompass a wide 
range of theories including Newtonian gravity and M0Ne13. 
We therefore, parametrize the gravitational force (F) at a 
radius r as, 



F oc 



(1) 



where 7 is a free parameter which is to be determined from 
observations. Here, 7 = 2 gives back the simple Newtonian 
case, while 7 = 1 gives the familiar MOND case in the low 
acceleration regime relevant at galactic scales. 

To relate the gravitational forc e to a rotation veloc- 
ity, on e needs a prescription for Md . [Garcia- Appadoo et al] 
(|2009l ) points out that the observed dynamical mass scales 
as cube of the optical radius. [Lpeb & Weincr (2011) have 
shown that dark matter particles interacting through a 
Yukawa potential could provide a natural explanation for 
a characteristic densit y in dark matter dominated h alos. We 
exploit this observed I Garcia- Appadoo et al.[[2009l ) relation 
Afd oc _R^, where R is the size of the Hi disk, to express the 
rotation velocity as. 



V (xR 



(2+7)/4 



(2) 



20031) reported that Mhi a nd surface areas of ADBS 



Rosenberg fc Schneidei] [200d ) galaxies are consistent with 



Thus the dependence of the asymptotic rotation velocity on 
disk size is a function of 7 and therefore is a probe of the 
effe ctive gravitational f orce law in the low acceleration case. 

[Chakrabortil (|201l[ ) has argued that the observed corre- 
lation between Hi mass and surface area may be the result 
of self regulated star formation. Even when a disk galaxy is 
driven away from this fundamental line, sel f regulation of the 
poros ity (driven by supernova feedback jMcCrav fc Snow 
'J97SJ) of the ISM and gravitational instability ( Safronov 
il96dt [Toomr3[l964) of the star forming disk, can control 
the star formation rate. This ensures that the star forming 
disk is alw a ys m arginally unstable. Under these consider- 
ations [MS (|1997[ ) has shown that the star formation rate 
(SFR) scales as, SFR oc V^^^, where V is the asymptotic 
rotation velocity of the disk. 



^ MOND can be thought of as a modification of either the inertial 
term (_F = fima) or of the effective gravitational force term (_F = 
^Newton/M) by fi{a/ao), where ao is the acceleration scale below 
which the dynamics is modified. At the low acceleration regime, 
relevant for disk galaxies, th ey result in the same rotation curves 
dSanders fc McGauBhll2002l'l . 
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Figure 1. Mhi vs A data fo r ADBS galaxies (cyan balls) from 
iRoaenberg fc SchneideJ l|2003') . Statistical uncertainties are com- 
parable to the sizes of the balls. Blue solid (and dotted) line shows 
the Newtonian+KS model (and 2a scatter region) with 7 = 2 and 
/3 = 1.4 put in Eqn. E] Red dashed line shows the MOND+KS 
model with 7 = 1 and /3 = 1.4. 

We can now express the surface density of the SFR 
(Esfr) as EsFR = SFR/(7r7i^) oc V^I'^Bt'^. Substituting 
for the rotation velocity (V) as a function of R from Eqn.[2l 
we have 

ESPR oc ii{5(2+7)/8)-2 ^ ^(5^-6)/8^ (3) 

This relates the surface density of SFR to the size of the 
galaxy. However, the Esfr is also related t o the surface den - 
sity of cold gas according to the KS law ijKennicuttll 19981 ). 
If the cold gas in a galaxy is dominated by Hi, then we 
can use the KS law to express the SFR surface density as 
SsFR oc S^j, where iKennicutti (1998) gives the power law 
index fi a value of 1.4 ± 0.15. 

Eliminating Esfr, we get the Hi surface density as 
E^j oc , Integrating the surface density over the 

size of the galaxy, we get 

Mhi = ^J?'Ehi oc 7?((5^-6)/{b«)+2. (4) 

Expressing this as a function of the surface area, we get 

Mhi = Cav1«^^~'''/<^''''»+\ (5) 

where Ca is a constant of proportionality. This general- 
izes the result (Mhi oc A^'^^ for Newtonia n gravity and KS 
law slope of 1.4) from IChakrabortil (j201ll ) to arbitrary KS 
law slopes and generalized effect ive gravitational force laws. 
iKhedekar fc Chakrabortil |2oi3) have shown that such a re- 
lation depends only weakly on the metallicity of the galaxy. 
The power law relates the Hi mass of a galaxy to its size 
and can be directly compared to the data (See Fig. [l|. 

We can now use the scaling in Eqn. [2] between V and R 
to find Mhi as a function of the galaxy size. This gives us 
the Hi mass as 

Mhi = CvI^(^^+i^^-^'/(^''(^+^», (6) 

where Cv is another constant of proportionality. This shows 
that the observed slope of the HITF relation is a probe of 
both the KS law and the law of gravity. For Newtonian grav- 
ity and KS law slope of 1.4 we have Mhi oc y'^-^^ , This can 
now be compared directly with the data (See Fig. [2]), but a 



Figure 2. A/hi vs inclination corrected W20 data for galaxies 
from Arecibo (magenta balls) and GMRT (cyan balls). See text 
for sources of data. Statistical uncertainties are comparable to 
the sizes of the balls. Blue solid (and dotted) line shows the 
Newtonian+KS model (and 2u scatter region) with 7 = 2 and 
/3 = 1.4 put in Eqn. [6] Red dashed line shows the MOND+KS 
model with 7 = 1 and /3 = 1.4. 

prior knowledge of the KS law slope will be required to de- 
termine the effective gravitational force law. This shortcom- 
ing can be overcome by simultaneously comparing Eqn.[6]to 
the HITF data and Eqn.[5]to the Hi mass vs area data. The 
ratio of the two slopes is a function of 7, but independent 
of fi. Hence, determining both slopes from observations can 
break the y — P degeneracy and determine 7 irrespective of 
our prior knowledge of /3. Being able to explain both ob- 
served relations for plausible values of these parameters is a 
crucial reality check for our model. 



4 CATALOGUE OF GALAXIES 

The Arecibo Dual-Beam Survey l|Rosenberg fc Schneider! 
boOOi' l (ADBS) has conducted a "blind" survey of ~ 430 
deg^ of sky and detected the Hi signal in 265 galaxies. 
Interferometric mapping of 84 galaxies was carried out 
with the NRAO's Very Large Array (VLA) as most of the 
ADBS galaxies were unresolved at the resolution of the 
Arecibo. Accurate sizes of 50 of them were compiled by 
[Rosenberg fc Schneider] (|2003| ) to tabulate Mhi and A (Hi 
cross section with column density above 2 x 10'^'' cm"'^) of 
individual ADBS galaxies. The observed sizes span 3 or- 
ders of magnitude (0.5 > log(yl/kpc^) > 3.5) and provide 
the data for the Hi mass vs area relation. For the follow- 
ing analysis, we have used only objects with large projected 
areas (log(j4/kpc'^) > 2) as the smaller galaxies may have 
systematically uncertain sizes limited by the resolution of 
the present surveys. Fig. [T] presents the observed data and 
com pares it with th e fiduci al models. 

(Curovich erahl (|2010l ) presented a sample of local 
{D < 60 Mpc) field galaxies with accurate 21 cm obser- 
vations to c onstruct a BTF r elation. We s elect the galax- 
ies listed bv lSakai et al] (|2000l ) and used bv lGurovich et all 
(201 01). All these galaxies have accurate Cepheid distances 
(Sak ai et al1l200Gl ) determined using the Hubble Space Tele- 
scope Ho Key Project and 21 cm Arecibo observations from 
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Figure 3. The hotter (cooler) colours denote regions with higher 
(lower) mean likelihoods. The solid (dotted) lines denote the 
la-/2a- confidence limits of the marginalized (mean) likelihoods. 
The figures on the right include a Gaussian prior on /3, while the 
figures on the left do not include such an external prior. Note that 
the value of 7=1, corresponding to the case of the MOND gravity, 
is completely ruled out, while the Newtonian gravity (7=2) gives 
an excellent fit to the data. 



iGiovanelli et al.1 (|l99it ). For all these galaxies we use Mm 
and W20 (the width at which th e Hi profile drop s to 20 % of 
its maximum value) as listed bv lGurovich et al. (|2010|). in- 
clinat ion correcte d following the method of Giovanelli et al.l 
(|l997l ). Following ITuIIv fc FisheH l|l977l ) this is commonly 
used as a proxy for rotation velocity in radio astronomy q 
This provides the data for the high mass range in the HITF 
relation. 

To accurately determine the slope of the HITF relation, 
this catalogue needs to be supplemented with data on low 
mass galaxies as well. The Fa int Irregular Galaxies GMRT 
Survey (FIGGS) conducted bv lBegum et al.l l|2008bl ) aims at 
characterizing the neutral ISM properties of faint, gas-rich 
dwarf galaxies through 21 cm GMRT observations. This pro- 
vides us with t he Mm and W20 (aga in inclination corrected 
following Ref. (|Giovanelli et al.lll997l )) measurements at the 
low mass end of the HITF relation. At the low velocity end, 
line broadening due to velocity dispersion from turbulent 
motion in the Hi, may play a significant rol e. We correct for 
this e ffect following the prescription from iTuUy fc Fouaud 
1 19851 ) B Fig [2] displays the observational HITF data and 
compares it with the fiducial models. 



^ It is possible to infer a rotation velocity Viot by fitting, say a 
tilted ring model (with a poorly determin ed inclination ) to the 
Hi data cube, but this is model dependent. iBegum et al~ 1 l l2008al ') 
confirm that W20 correlates better with mass indicators than VJ-ot . 
In the rest of the analysis we replace Vrot with its observational 
proxy W20- 

^ This may intr oduce a syst ematic uncertainty. However, accord- 
ing to Eqn. 5 of lSilkl |200l), a 10% variation in the gas density, 
leads to a ~ 1% change in the gas dispersion velocity. This will 
introduce a < 2% systematic error in the corrected velocity width, 
which is small compared to the statistical uncertainty of ~ 10%. 



5 MCMC ESTIMATION 

In order to obtain an effective law of gravity at the galactic 
scales using the data on Mm, A and v, we consider the two 
scaling relations - Mm vs A and Mm vs v, as expressed 
through Eqn. [5] and [6] respectively. We compute for the 
the data as a sum, — Xa + xt, one for each for the 
relations - Mm vs A and Mm vs i;, respectively. For each x^- 
The weight, a, is computed by adding the intrinsic scatter 
and observed error in quadrature. 

We then compute the likelihoods, exp(— in the 
four dimensional parameter space of {Ca, Cv, P, 7}. We 
perform Markov Chain Monte Carlo (MCMC) runfl in the 
using these likelihoods. Of these parameters, we eventually 
marginalize over the nuisance parameters Ca and Cv- Note 
that we do not assume a completely deterministic form for 
the KS law, but require only a simple power law relation, 
which is eventually marginalized when quoting our final re- 
sults on 7. 

The results of our likelihood analysis is shown in Fig. 
[3l The top figures displays the la and 2a constraints in the 
/3 — 7 plane. The top left figure shows the likelihoods without 
any external prior on /3, while the top right figure shows 
the results after inc lusion of a Gauss ian prior of /3 = 1.4 ± 
0.15 as obtained by iKennicuttI (|l998l ) through independent 
observations. The corresponding figures at the bottom show 
the marginalized (over /3) 1-D distribution on 7. 

Our MCMC exploration of the parameter space, with- 
out the priors on the star formation law, reveals that even 



though P = 2.01 



+4.23 



is not well constrained, the effective 



law of gravity is highly constrained with 7=2.06±0.34, con- 
sistent with Newtonian gravity (7 — 2) but clearly ruling 
out MOND (7 = 1). At the same time, the best fit region 
in the /3 — 7 plane also i ncludes the regio n with KS law 
slope of /3=1.4±0.15 from iKennicut^ j 19981 ). Hence our re- 
sult confirms that the observed HI scaling relation is not 
only consistent with Newtonian gravity but also with the 
KS law prescription for star formation q We find that both 
with and without an external prior on /3, our model can 
conclusively rule out the MOND model of gravity; the con- 
straints on 7 get slightly better after including the prior on 
p. This shows that our results are insensitive to the prior 
used on the KS law, as expected from Section O The likeli- 
hood analysis yields values of 7=2.08 ± 0.29 after including 
an external prior on /3. 



6 DISCUSSIONS 

TF relations are being increasingly discussed (and also crit- 
icised) in the context of weak field tests of gravity. In this 
letter we have suggested the simultaneous use of two Hi 



* An MCMC algorithm is useful for sampling the likelihoods from 
a high dimensional probability distribution based on constructing 
Markov chains that have the desired di stribution as their equ ilib- 
rium distribution, i.e. exp(— x^/2); see lLewis &: Bridlj l l2002t) for 
more details. 

^ iRovchowdhurv et al. 1 1I2OIII ) have found that when comparing 
Sfii to the Ha emission, the KS law is steeper (/3 ~ 2), which 
is also consistent with our result. A summary of the impact of 
different star form ation laws on galaxy formation is given by 
iLagos et all | |2010|'I . 
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scaling relations for this purpose. Our test is free from some 
of the systematics and selection bias that affect the exist- 
ing methods. Our model incorporates, for the first time, the 
physics of supernova feedback driven hot gas into this anal- 
ysis. We have shown that a simple model of self regulated 
star formation accounts for the observed Ifl scaling relations 
in disk galaxies. 

The crucial ingredients which make the model agree 
with the data are the star formation law and the effective 
law of gravity in the weak field limit. Earlier work shows 
that rotation curves and single scaling relationships, such as 
the B TF, can be e xplained both w ithin the DM (jLarsonI 
l20ld) and MOND (|McGaughl boilT ) paradigms. Hence, a 
new method was required to use galactic dynamics for a 
weak field test of gravity. Our work shows that the HITF 
and the M-vs-A relations taken together can constrain the 
effective law of gravity relevant on the galactic scale. The 
data rules out MOND at more than 99.9% confidence at 
galactic scales within the scope of this model. We show that 
our model suffices to explain the data using a simple KS law 
prescription for star formation and Newtonian gravity. This 
is a triumph for the DM paradigm with Newtonian gravity 
in the low acceleration regime. 

With Extended-V LA follow up of ALFALFA 
jGiovanelli et al.l l2005f ) detected galaxies and the the 
advent of high sensitivity, high resolution, next generation 
radio telescopes such as the Square Kilometer Array, the 
number of Hi detected galaxi es is set to grow rapidly 
IIKhedekar fc Chakrabortil I2OIII ). Our study implies that 
the soon to be observed scaling relations between Hi masses 
of these galaxies, their rotation velocities and sizes will 
provide an excellent probe for the nature of gravity and 
dynamics in the weak field limit. 
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